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Abstract 
Q_[ 

D ■ The first two non-trivial moments of the distribution of the topological charge (or gluonic 

C/3 . winding number), i.e., the topological susceptibility and the fourth cumulant, can be com- 

puted in lattice QCD simulations and exploited to constrain the pattern of chiral symmetry 
breaking. We compute these two topological observables at next-to-leading order in three- 
flavour Chiral Perturbation Theory, and we discuss the role played by the r\ propagation 
in these expressions. For hierarchies of light-quark masses close to the physical situation, 
we show that the fourth cumulant has a much better sensitivity than the topological sus- 
ceptibility to the three-flavour quark condensate, and thus constitutes a relevant tool to 
determine the pattern of chiral symmetry breaking in the limit of three massless flavours. 
We provide the complete formulae for the two topological observables in the isospin limit, 
and predict their values in the particular setting of the recent analysis of the RBC/UKQCD 
collaboration. We show that a combination of the topological susceptibility and the fourth 
cumulant is able to pin down the three flavour condensate in a particularly clean way in 
the case of three degenerate quarks. 
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1 Introduction 



An accurate description of low-energy QCD remains a very challenging task for the theorists. One of 
its prominent features, the spontaneous breakdown of chiral symmetry leading to Goldstone bosons 
identified with the pseudoscalar ir, K, rj mesons, can be described through an effective theory, Chiral 
Perturbation Theory (xPT) [HE]- However, this tool based on the symmetries of the underlying 
theory does not allow one to determine the precise pattern of chiral symmetry breaking. A very 
helpful complement to the effective field theory approach comes from numerical simulations of Quan- 
tum Chromodynamics on the lattice, as they provide a way to determine the low-energy constants 
(LECs) encoding the pattern of chiral symmetry breaking through the momentum and quark-mass 
dependence of various observables related to the dynamics of pseudoscalar mesons [3] (and conversely, 
chiral expansions can be exploited by lattice collaborations to extrapolate their data down to physical 
values of the light quark masses) . 

Another interesting way to determine non-perturbative features of the QCD vacuum consists 
in exploiting the connection between the fluctuation of the winding number, or topological charge, 
defined from the gluonic strength tensor G as 



and the breaking of the anomalous axial rotation £7,4.(1), related to the dynamics of the r\ and rf 
mesons [IH7]. Indeed, the fluctuations of the topological charge correspond to long-distance dynamics 
encoded in \PT: and can be monitored during the Monte Carlo evaluation of the path integral in 
lattice simulations. The study of topological observables is thus able to provide additional information 
on the pattern of chiral symmetry breaking (see ref. [8] and references therein). 

In a previous article [9J, we have discussed the potentiality of the variance of the topological 
charge, also called topological susceptibility and measured in lattice simulations, to determine the 
pattern of chiral symmetry breaking. We have used data available from the RBC and UKQCD 
collaborations (in combination with information on the spectrum of light pseudoscalar mesons) in 
order to constrain this pattern. We have shown that for simulations close to the physical point, the 
topological susceptibility only probes the quark condensate in the chiral limit of two massless flavours 
(m u ,md — > 0, m s kept at its physical value), but not the condensate in the three-flavour chiral limit 
{m u ,md,m s — > 0). Indeed, in ref. [9], the pattern of three-flavour chiral symmetry breaking was 
mainly constrained by the spectrum of light pseudoscalar mesons, while the topological susceptibility 
played only a minor part in the analysis. The two-flavour condensate has been measured through irn 
rescattering in various low-energy processes, such as [10H16] and K — > 3tt decays |17H21j . But 
this condensate can be significantly higher than the three-flavour quark condensate in the presence 
of large contributions from sea-quark ss pairs to the chiral structure of QCD vacuum [22] , a scenario 
supported by several phenomenological analyses (irK scattering [23J12I], scalar form factors [2"5Tf2"5] , 
lattice data [9,29j). 

In the present work, we will investigate a slightly more sophisticated topological observable, 
namely the fourth cumulant C4 of the gluonic winding number. This quantity can be measured on 
the lattice [8l l30ff3l] . and we will see that it does not suffer from the same shortcomings as the 
topological susceptibility to determine the three-flavour chiral condensate for simulations close to 
the physical point. In a similar way to the topological susceptibility [9,35,36j, the fourth cumulant 
C4 can be computed at NLO in xPT, and this will be the focus of the present article. We will 
compute the fourth cumulant C4 in two different ways. In Sec. El we will determine it by considering 
the effective potential for constant sources and obtaining its dependence on the vacuum angle 9, 
similarly to Ref. [35]. In Sec. El we will express C4 as the 4-point correlator of the winding number 
density uj at vanishing incoming momenta, and analyse the correlator at one loop in \PT using 
Feynman diagrams, so that we can single out the role of the 77 propagation. In Sec. [H we analyse 
briefly the potentiality of the fourth cumulant to extract the pattern of three-flavour chiral symmetry 
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breaking from lattice simulations close to the physical point, with a much greater sensitivity than the 
topological susceptibility. We will also focus on the expansion in two-flavour xPT, before drawing a 
few conclusions. Two appendices describe details of the computation using Feynman diagrams, as 
well as the reexpression of NLO low-energy constants in terms of physical quantities. 



2 Derivation through the effective potential 
2.1 Distribution of the topological charge 

Let us first define these two quantities related to topological aspects of QCD [6,8,37J. The partition 
function of QCD in a ^-vacuum with a given quark-mass matrix M is defined as 

z(M,e) = J2^ iQ6 z Q (M), (2) 

Q 

where Zq is the partition function for a fixed topological sector, with a fixed winding number (or 
topological charge) Q, which is an integer in presence of fermion fields in the fundamental represen- 
tation. If we work in a finite volume V with an Euclidean metric, we can define the vacuum energy 
density as 

e vac (M,9) = - ±- log Z(M, 9), (3) 



V 



and its second and fourth cumulants 

5 2 e vac (M,9) 



X 



56 2 



5 4 e vac (M,9) 

c 4 = 



6»=0 



89 4 



(4) 

0=0 



Both quantities can be interpreted as quantities describing the distribution of the topological charge 

Q 

(Q 2 )e=o [(Q 4 )-3(Q 2 )Vo 

x = — — C4 = v • (5) 

X, related to the variance of the distribution, is the topological susceptibility and has been studied 
in detail in Refs. [6l|9l[35j[36]- C4 measures the kurtosis of the distribution, i.e., its more or less 
peaked nature (a Gaussian distribution has a vanishing kurtosis). Both quantities can be defined non- 
ambiguously in terms of Green functions of scalar and pseudoscalar densities at vanishing momentum, 
so that they are renormalisation-group invariant [38J. 

For large volume and small quark masses, the partition function is dominated by the lightest 
states of the theory, i.e., the pseudoscalar Goldstone bosons, and can thus be analysed in terms of 
Chiral Perturbation Theory. As we deal only with constant source terms (quark masses and vacuum 
angle), we can focus on modes with vanishing momentum, i.e. the matrix U collecting Goldstone 
bosons become independent of x 

Z{M, 6) = J [dU] exp[-VC(U, M, 9)] = exp[-VW(M, 9)] , (6) 

where C is the xPT Lagrangian (which can be considered for U constant) and W is the (Euclidean) 
effective potential. As shown in Ref. [39], the effective theory of QCD in a finite volume with periodic 
boundary conditions amounts to the same Lagrangian as in the infinite volume. The properties of 
this partition function have been extensively studied in Refs. [6 ll40j . and in particular the distribution 
of the winding number Q according to the leading-order (LO) chiral Lagrangian. 
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2.2 Structure of the one-loop generating functional 

We can follow the arguments of Ref. to derive the expression of C4 in an elegant way using the 
effective potential of the theory Eq. ([6]). The latter can be obtained from the founding work of xPT, 
ref [2], where one can find the one-loop generating functional Z up to 0($ 4 ) in the Minkowskian 
metric, with the following counting of the sources a M ~ p = 0(<3?) and ~ s — M = 0(Q 2 ), but in 
the absence of a source term for the vacuum angle. This is however sufficient for our purposes, as a 
constant source term 9 can be introduced in the effective potential via an anomalous Ua(1) rotation 
which leaves invariant the generating functional 

Z(s + ip = s + ipo, 6 = e ) = Z((s + ip )e ie °/ N , 0) , (7) 

as discussed in detail in ref. [2J. The generating functional has been extended later for an arbitrary 
number of light flavours N in refs. |4HI42j. and will be used here to determine the contribution of 
order O(0 4 ) in the effective potential. We set the sources 

= a M = 0, (s + ip) = Me ie ' N , M = diag(mi ,...m N ). (8) 

It is easy to determine the classical solution to the LO equation of motion under the form 

U = (e^...e^), 5>i = 0, (9) 

j 

which should minimise the LO Lagrangian 

£1 = FjjB N ^ m j cos(0j) , (10) 
j 

with 4>j = 9/N — ctj with the condition ^ ■ cf>j = 6 to be fulfilled 0. This Lagrangian involves the 
pseudoscalar decay constant Fjy and the chiral condensate S(A r ) = —F^B^ in the chiral limit of N 
massless flavours Q 

Using a Lagrange parameter, we can determine the minimum as an expansion in powers of 9 



mi 



m \ I m \ m 3 



j+G(9 5 ), (12) 



rrii J \rrii J mi 3 ! 

where we have introduced the notation for the various sums and harmonic means 

™W = — L_ } s lk] = J2m^, fh = m [1] , s = s [1] . (13) 



E 



J m 



We can then plug this expression into the one-loop generating functional Z evaluated at the point 
Eq. ([8]). The partition function gets contributions from two different terms Zt and Z u , collecting 



1 For simplicity, we have chosen to apply the same phase e l9 ^ N on each of the quark masses in Eq. (|5). but we could 
have chosen any redefinition of the form 

«m = dp = 0, (s + ip) = diag(mie™ lS , . . .mjve™" 9 ) , Wj = 1 . (11) 

i 

The rest of the demonstration is unchanged, provided that one adapts the definition of <j>i — Wi9 — a.j. In particular, 
one notices that the effective potential is only defined in terms of cj>i (and not 6 and Oi separately), so that our results 
for the topological susceptibility and the fourth cumulant are independent of the choice of weights lOj, as expected. 

2 Fjv and Bn would be denoted Fo and Bo in the case of N — 3 [2], and F and B in the case of = 2 x?T 
respectively [l]. 
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tadpoles and unitarity contributions respectively [2,42j (there is no contribution from the anomalous 
part Za) 



z t /v 



z u /v 



with 



a 



PQ 



u 



N 



2(N 2 - 1 



-F 



N 



o 2 

M P 
32tt 2 



log 



o 2 ' 

M P 
M 2 



Opp + \§B\L\. N {p) y^mj cos i 



-l§B 2 N U 7 . N (n) I my sin <ftj j + 8BjfLl. N (fi) y^ m| cos 20 
lJ2 I dx J PQ( x ) a PQ a . 



^QP 



P,Q ' 



(e 



ioi/2 



Ja n /2 



) , A = (mi (cos 0i — 1), . . . mAr(cos 0tv — 1)) 



(14) 



(15) 



^<{T P ,Tt }(„Mue-^ + - 5 PQ M P = ?f({T P , T&A) , (16) 



(17) 



The diagonalisation of the mass term in the LO equation of motion defines Mp, the leading-order 
contribution to the masses of the N 2 — 1 pseudoscalar mesons and Tp, the linear combinations 
of the generators of the SU(N) flavour group (corresponding to Gell-Mann matrices for N = 3). 
J r corresponds to the one-loop scalar integral with two mesons P and Q, and Lq,Lj,Ls denote 
Low-Energy Constants (LECs) renormalised at the scale //, entering the NLO chiral Lagrangian Li 
expressed in the basis of operators detailed in Ref. [2] for N = 3 xPT (the same basis holds for an 
arbitrary number of flavours, up to the introduction of an additional LEC Lq that does not enter the 
present computation |41U42j ). Using Eq. (|12|) . we have the expansion of the trigonometric functions 



,-• cos ' 



j 
j 

y^ m 2 cos 2cj)j 



sin * 



m T + ^W2l + 0( * 



Nm6 - N 



m 



4 /)3 



[3] 6 



m 



+ 0(9 6 ), 



- 2Nm 2 2 + In^-t9 a + 0{6 & ) , 

3 TTiL^J 



(19) 
(20) 
(21) 



as well as the expansion of a x 



mi 



mi) 2 I rm 



m 



in 



and the scalar integral at vanishing transfer momentum 



/ 



o 2 

M f 



31- 

mi 



o 2 



dx J 



-2k 



PQ 



1 M P log ^-Mq log 
16^2 2 72 

M P - Mq 



Mr 



(22) 



(23) 



so that all the elements in Z% and Z u are simple functions of the quark masses and the chiral LECs. 



3 For N = 3 xPT, one has 



Mn= 2mB , 



M K = {m + m s )B , 



o 2 2 

M v = n( 2m = +m)B . 



(18) 
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2.3 NLO expression of the topological susceptibility and the fourth cumulant 



We are now in a position to determine the one-loop expression of the first two cumulants of the 
winding number. At the order 0(9 2 ), needed for the topological susceptibility %, on has the following 
contribution 

2 



Zt/V - 



° 2 ■ ., 

B N F 2 N m + B N ^^^g- 

PA 



{T P ,T],}i 



(24) 



-32B 2 N L r 6 . N (n)sm - 32B 2 N L r 7 . N {^)N z m I - 32B z N L r 8 . N (fi)N 



2-2 



?2 rr 



m 



ZJV -> 0, 
Using the summation formula 



4(iV 2 



N 



J%3 > 



(25) 



(26) 



the expression of the topological susceptibility at one loop given in Refs. [9, 35, 36J is recovered in a 
straightforward way 



X = B N F 2 N m-B N J2 ^ log 

PA 



M P m z 
fj, 2 nii 



{rp,4} ii +32B^LS ;JV (/i)sm+325^iV[JVL?. JV (/i)+LS. JV (/i)] 



m 



(27) 

Let us notice that Ref. [35] obtained this result by determining the classical solution Eq. (fT2j) cor- 
responding to the minimum of the chiral Lagrangian up to next-to-leading order C\ + £2- This is 
actually an unnecessary complication, since the expression of the one-loop effective potential given 
in Ref. [2] is precisely designed to require the classical solution from the leading-order Lagrangian 
only. Indeed we recover the expression NLO of the topological susceptibility without being forced to 
perform the elaborate minimisation of Ref. |35j . 

In the isospin limit m = m u = for N = 3 flavours, one has 0E51E6] 

nopole B F 2 m m s , 32mm s B 2 Ll(fi){2m + m s ) t 96m 2 m 2 B 2 [3L 7 + L§(//)] 



X 



+ v ~ + 

m + 2m s m + 2m s 

o 2 

M n mm s B 2 (m + m s ) 2 
2 



3B 2 m 2 m 2 

— log 

87r 2 (m + 2m s ) 2 \i 



8ir 2 (m + 2m s ) 2 



log 



[m + 2m s ) z 

o 2 

M K 



(28) 



o 2 



mm s B${2m + m s ) pol pole 

727r 2 (m + 2m s ) S /i 2 ^ X * 



where c£J° pole is a remainder collecting higher-order (HO) contributions, which starts at 0(m 2 ). 



At the order <3(6> 4 ), the one-loop generating functional reads 



1 

24 



m 



[3] 



m 



+32B^LU(/x) 



o 2 

64vr 2 



log 



o 2 
tfp 

M 2 



{Tp,^}^ — 



m 



[3] 



m 



. rrii 



(29) 



??? 



[3] 



+ 128B 2 N L r 8 . N (fi)N 



m 



[3] 



ZJV 



6 4 3Bjj _ 4 

— m 

24 4 



(30) 



4 The explanation of the superscript "no pole" is given in Ref. [5] and will become clear in the following. 
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which yields the one-loop expression of the fourth cumulant 

o 2 2 



- 3 / - \ 2 

4^-3 — 



4 \ - 5 



(31) 



-32^.^) + 3m 2 ) - 128B 2 N N[NL r 7] M + ^ ;JV ( M )]^ 

PQ \ i 7 \ i J 



. 3-B^r _ 4 

H -r-m 



To our knowledge, Eq. (|3ip is the first NLO computation of the fourth cumulant in xPT. Let us 
add that this formula agrees with the leading-order result presented in Ref. [35]. Eqs. (f27"|) and (f3T]) 
could be used in principle to derive the expression of x an d C4 for an arbitrary number of flavours 
N. However, one should emphasise that the notation and definitions of the LECs in these equations 
are derived from the basis used in the three-flavour case that will be our main focus. In Sec. we 
will discuss how these formulae must be rewritten to match the usual definition of the LECs in the 
case oi N = 2 x?T built around the two-flavour chiral limit m u = m^ = 0, involving only soft pions 
as dynamical degrees of freedom pQ. 

In the case of three flavours, it is straightforward to check that our NLO expression for C4 is 
indeed scale independent, even in the presence of strong isospin. In the isospin limit m = m u = m^, 
one can find easily the expression for C4 

no pole _ B FQmm s (m 3 + 2m 3 ) 64£? 2 mm s Lg(/i) (m 4 + 2m 3 m s + Qm 2 m? s + 8mm 3 s + mf) 
° 4 ~ (m + 2m s ) 4 (m + 2m s ) 4 ^ ' 

o 2 

384S 2 m 2 m 2 (3L 7 + Lg(/x)) (m 3 + 2m 3 ) 3m 2 m 2 5 2 (m 3 + 2m 3 ) J M v 
(m + 2m s ) 5 27r 2 (m + 2m s ) 5 ^ /x 2 

o 2 

mm s BQ (m 3 + 9mm 2 + 2m 3 ) (m + m s ) 2 Mjf 
87r 2 (m + 2m s ) 5 ^ /i 2 

o 2 

BQmm s (m 4 + 2m 3 m s + 6m 2 m 2 + 8mm 3 + m 4 ) M„ S 2 m 2 m 2 (l3m 2 + 22mm s + 37m 2 ) 
+ 36vr 2 (m + 2m s ) 4 log 7^" + 24vr 2 (m + 2m s ) 4 

. no pole mo pole 
"r c 4 a c 4 J 

where dc° P ° le is a remainder collecting HO contributions, starting at 0(m 2 ). 

The above formula features chiral logarithms of different origins. Some of them come from tadpole 
contributions (Zt), whereas others stem from the loop function Jpg taken at vanishing momentum 
transfer as indicated in Eq. (j23|) (Z u ). At this order, one could in principle redefine the argument of 
(some of) the logarithms in order to have physical masses instead of leading order masses. This change 
would not affect the remainder of the NLO expansion, apart from a redefinition of HO remainders. 
Following the discussion in Refs. |24 l l29 l l4"3]). we will consider either the above prescription where 
all the logarithms have LO masses as their arguments, or the one where we take physical masses for 
the unitary logarithms but LO ones for the tadpole logarithms El If we perform this separation, we 



5 This choice is made to meet two different criteria. First, we want to avoid redefinitions which are not motivated 
by physical arguments, which invites us to keep LO masses as the argument of the tadpole logarithms, as we have no 
intuition on how these logarithms will be affected by higher and higher-order contributions in the chiral expansion. 
Second, as illustrated in Refs. [291143] . the unitarity logarithms Eq. Q23p arise in the momentum dependence of form 
factors and scattering amplitudes of pseudoscalar mesons and will induce unphysical divergences in the case of a 

o 2 

vanishing quark condensate where Mp— > 0, supporting a reexpression of the arguments of the unitary logarithms in 
Eq. (|23p in terms of physical masses. 
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obtain 



c no pole = B Q F 2 mm s (m 3 + 2m 3 s ) _ 384B 2 m 2 m 2 [3L 7 + L§(//)] (m 3 + 2m 3 ) 
4 (m + 2m s ) 4 (m + 2m s ) 5 

64BqLq (fj,)mm s (m 4 + 2m 3 m s + 6m 2 m 2 + 8mm 3 + m 4 ) 
(m + 2m s ) 4 

o 2 

B^m 2 m 2 a {13m 2 + 22mm s + 37m;?) JVf^ 3B 2 m 2 m 2 (4m 3 - 3mm 2 . + 2m 3 ) 

+ 247r 2 (m + 2m s ) 4 + l0g 7^ 8vr 2 (m + 2m s ) 5 

o 2 

j Mx B 2 mm s (m + m s ) 2 (m 3 — 3m 2 m s + 3mm 2 + 2m 3 ) 



(33) 



/i 2 87r 2 (m + 2m 



o 2 

M~ B^mm s (2m 4 — 8m 3 m s + 13mm 3 + 2m 4 ) 
+ l0g 7^ 72vr 2 (m + 2m s ) 4 

M 2 9 J B 2 m 2 m 4 M| 3B^m 2 m 2 (m + m s 

+ 0g "^ 2 "8vr 2 (m + 2m s ) 4 + ° g 8vr 2 (m + 2m s ) 4 

, i „ M ^ ^gm 2 m 2 (2m + m,) 2 no pole no pole 
+ 1 ° g /x 2 24vr 2 (m + 2m s ) 4 + ° 4 ^ 



It is obvious that the HO remainder pole has absorbed the redefinition of the argument of the 
logarithms. In the following section, we will show the computation for C4 using Feynman diagrams and 
performing this distinction, keeping in mind that we can always revert to the prescription in Eq. (j32|) 
by setting LO masses in the argument of all the logarithms. For both prescriptions Eqs. (j32|) and 
(|33|) . we see that C4 vanishes in the limit where m = or m s = 0. This is expected since the effective 
potential becomes independent of 9 in the limit where at least one of the quark masses vanishes, as 
the vacuum angle can be rotated away through an anomalous Ua{1) rotation (see Refs. and 
references therein for a more detailed discussion). 



3 Analysis in terms of Feynman diagrams 
3.1 Combinatorics 

It is also possible to derive the value of C4 using the formalism of Feynman diagrams in N = 3 xPT. 
It is obviously completely equivalent to the previous approach in terms of the one-loop generating 
functional, but it allows one to separate the contributions coming from the propagation of mesons. 
Computing the fourth derivative of the generating functional with respect to 9, we obtain 



5 A Z 



59* 





J [dU]e iL 







5 4 L . 5 3 L 5L 



'59* 



59 3 59 



6 \59 2 



6i 



\59 2 



f8L 
\59 



+ 



SL 



,(34) 



J e=o 



with L = J C. At leading order, only the LO chiral Lagrangian L\ is needed, and more precisely, 
as all incoming momenta vanish, the mass term with Me t6 l N , see Eq. (|l(jp . The five tree diagrams 
generated are indicated in Fig. [lj in the same order as the derivatives in Eq. (|34p - the third term 
in Eq. (|34|) cannot yield any contribution to C4 at tree level, whereas the fifth term yields both 
diagrams 4 and 5. In Fig. [H the number inside each circle indicates the number of derivatives with 
respect to 9 applied to L\ to obtain the corresponding ^-induced vertex, whereas the four-leg central 
vertex in diagram 5 is obtained directly from C\ (taking the term proportional to -Bo, which has no 
derivatives) . 
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Figure 1: Tree diagrams contributing to C4. The 9 -induced vertices are denoted by circles containing the 
number of derivatives with respect to 9 applied to the mass term from the LO Lagrangian in order to obtain 
the corresponding vertex. 



It is easy to compute these diagrams, with the respective contributions in the isospin limit 

4 (m — m s ) 2 



jmo pole 



1 . , 8 (m — m s ) 2 

— (2m + m s ) H ^ -- 

81 V ; 81 m + 2m s 



27 m + 2m s 



(35) 



8 (m — 4m s )(m — m s ) 3 2 (m 



) (m + 8m s 



81 (m + 2m s ) 3 



81 81 (m + 2m s 



BoFQmm s (m 3 + 2m? s 
(m + 2m s ) 4 



+ 0{ml) 



adding up, as expected, to the LO result obtained in the previous section. In these diagrams, the 
meson lines are always coupled to at least one ^-induced source coming from a single derivative with 
respect to the vacuum angle (represented by a circle containing "1") 



Fo 2 B 



(M$) + 0($ 3 



(36) 



where U = exp(i<3?). In the isospin limit, one can see that the lines involving a single meson emitted 
from such a ^-induced vertex can correspond only to the propagation of an r\ meson (there would be 
an additional 7r° contribution out of the isospin limit). This explains the presence of denominators 

o 2 

with up to four powers of M v = 2(2m s + m)/3 in Eq. ([33]) . 

Moving to next-to-leading order, we can dress the tree diagrams by a) inserting a vertex from 
the NLO Lagrangian £2 inside one of the propagators, b) replacing a ^-induced vertex from L\ by its 
counterpart from £2, c) adding a loop to one of the ^-induced vertices from C\, d) adding a vertex 
from C\ to build a loop either as a rescattering diagram or a tadpole on one of the propagators. These 
operations are illustrated in the case of diagram 5 in Fig. [21 and the equivalent one-loop diagrams for 
diagrams 1-4 can be obtained straightforwardly. The number inside each circle (respectively square) 
indicates the number of derivatives with respect to 9 applied to C\ (respectively £2) to obtain the 
corresponding ^-induced vertex. In addition to dressing the diagrams in Fig. [lj one can also take 
more legs out of the ^-induced vertices from L\ and draw new diagrams, given in Fig. [3l Simple parity 
arguments show that ^-induced vertices containing an even (odd) number must have an even (odd) 
number of propagators attached to them. We provide the contributions of the different diagrams in 
App.[A] and one can check explicitly that the sum of all these diagrams yields Eq. (|33|) as expected. 



3.2 Isolating the rj propagator 

According to the previous analysis, the powers of (m+2m s ) in the denominator of Eq. (|33|) are related 
to the propagation of rj mesons with vanishing momentum, which are the only states coupling to the 
^-induced vertices in the isospin limit. Since we used NLO xPT, it is normal that the mass involved 

o 2 

in the r/-propagator is the LO value M v = 2(2m s + m)/3. For instance, the LO expression for C4 
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Figure 2: Scattering diagrams contributing to C4 at one loop and obtained by dressing the tree diagram 5, 
either by adding a tadpole loop (first two lines) or by replacing a LO vertex by its NLO counterpart (denoted 
by a square). The 9-induced vertices are denoted by circles (squares respectively) containing the number of 
derivatives with respect to 9 applied to the mass term from the LO (NLO respectively) Lagrangian in order to 
obtain the corresponding vertex. 



Eq. ([35]) should be understood as: 



~B F^2m 

?4 772 



ms) + 



16-Bq-Fq(to — m s ) 2 16-Bq(?71 — m s ) 2 (m + 2m s ) 



243M2 



243M4 



(37) 



+ 



64BqFq (m — 4m s )(m — m s ) c 
2187M6 



32^F 2 (m 



m s ) (m + 8m s ) 



6561M8 



+ 0(ml) 



However, if we go to higher and higher orders in Chiral Perturbation Theory, there would be tadpole 
and counterterm contributions to the propagator which would shift the propagator masses from LO 
to physical value. This is already the case with some of the contributions in diagrams 1-5, which 
arise when one inserts a tadpole or a NLO LEC in one propagator. These "double-propagator" 
contributions should actually be reabsorbed inside the mass of the 77 propagator to shift its pole from 

o 2 

M v to M 2 . 

This problem is all the more acute if we do not assume that three-flavour chiral symmetry 
breaking is triggered by large values of the chiral condensate and the pseudoscalar decay constant, 
but exhibits a more complicated (and interesting) pattern. In this case, the LO contribution to 
the 7] mass might be significantly different from its physical value, due to large contributions from 



9 






10 



Figure 3: Additional diagrams contributing to C4 at one loop and not derived by dressing tree diagrams. 

o 2 

NLO terms in its chiral expansion. In view of the large powers of jVf or M 2 involved (up to the 
fourth order for diagram 5), this distinction might be quite important, and suggests that one should 
identify the contributions due to 77 propagation in the previous computation and replace the LO rj 
mass by its physical value. This point of view is in agreement with the philosophy of Resummed 
Chiral Perturbation Theory [9,22,24,27-29,43-45], which is built to accommodate such patterns of 
chiral symmetry breaking and where physically-motivated redefinitions of the chiral expansions are 
performed in order to limit the size of higher-order remainders. 

In the present case, it is quite straightforward for each diagram to identify the NLO contributions 
that should be absorbed in the shift since the structure of the rj mass up to NLO is well known [2]. 
This yields the following expressions for the different diagrams, which should be summed together 
to obtain the NLO expression of C4 with a physical 77 propagator. 

• Diagrams 1 

1 B F 2 (2m+ms) + S ° m2 In ^ I B 2 (m + m s ) 2 M\ g| (m + 2m s ) 2 M 2 
--B F (2m + ms) + ^^ln- T + — 2 In + In (38) 

-^B 2 (m s + 2m) 2 L^) - ^-B 2 (m s + 2m) 2 L 7 - ^B 2 (2m 2 + m 2 ) 

• Diagrams 2 

o o 

16B$F$(m-m s ) 2 2B$m 2 (m - m s ) M 2 2B$ (m - m s ) (m + m s f M 2 K 

243M2 81vr 2 M2 ^ + 243vr 2 M2 /1 2 ( j 

o 

2Bq (m — m s ) (m — 4m s ) (m + 2m s ) M 2 
2187vr 2 M 2 n ~Jt? 

256i?g (m — m s ) 2 (2m + m s ) 256-Bq (m — m s ) 2 (m + 2m s ) 

243M| Li[ll) 729M| ^ 

2560^ (m - m s ) 2 (2m + m s ) 2560^ ( m - m s ) 2 (2m + m s ) 

+ 243M| Le(/i) + 243M| L7 

2560Sq (m — m s ) 2 (to + m s ) 
+ 243M| L * W 
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• Diagrams 3 



16-Bq (m — m s ) 2 (m + 2m s ) B^m? (m — m s ) (m + m s ) M 2 

243~mJ + 27vr 2 M 4 T 2 " { ' 

o o 

i?o (m - m s ) (m + m s f M 2 K (m - m s ) (m + 2m s ) (m 2 + mm s - 8m 2 s ) M 2 

162vr 2 M4 729vr 2 M 4 n 7?" 

512i?Q (m — m s ) 2 (2m + m s ) (m + 2m s ) 512Bq (m — m s ) 2 (m + 2m s ) 



243M4 ^ + 243M| ^ 

2040-Bq ("i — m s ) 2 (m 2 + mm s + m 2 ) r 2040-Bq ( m ~~ m s) 2 ("i 2 + mm s + m 2 ) 

L& ^ ] 8ll^ L? 

1024Bq (m — m s ) 2 (m 2 + mm s + 2m 2 ) 

-L 8 W 



81M4 



• Diagrams 4 



64SoF 2 (m - 4m s )(m - m s ) 3 AB^m 2 (Am - 13m s ) (m - m s ) 2 M 2 

2187M6 729vr 2 M6 n ( ^ 

O 

2B£ (m + m s ) (m - m s ) 2 (32m 2 - 99mm s - 23m 2 ) M% 

■In- ' 



10935vr 2 M6 ^ 2 

o 

ABk (m + m s ) (m — m s ) 2 (Am 2 — A\mm s + 64m 2 ) M, 

In 



19683vr 2 M6 " /i 2 

10245^ (m - m s ) 3 (m - 4m s ) (2m + m s ) 
729M| 

10245^ (m - m s ) 3 (m - 4m s ) (m + 2m s ) 

2187M| Ls( ^ 

40965^ (m - m s f (7m 2 - \lmm s - 14m 2 ) r 
+ 2187M| L%[jl) 

40965^ (m - mgf (7m 2 - llmm s - 14m 2 ) 
+ 2187M6 Lr 

20485^ (m - m s f (7m 2 - 9mm s - 28m 2 ) 
+ 2187M| Ls(M) 
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Diagrams 5 



o 

32B%F$(m-m s ) 4 (m + 8m s ) 2B%m 2 (m - m s ) 3 (5m + 31m s ) M 2 

6561M8 + 2187tt 2 M8 n 7^ ( ^ 

o 

2B% (m - m s ) 3 (m + m s ) (m 2 - 189mm s + 8m 2 ) ^ 
32805vr 2 M8 R 



2B% (m - m s ) 3 (m + 2m s ) (5m 2 + 47mm s - 160m 2 ) M 2 
+ 59049vr 2 M8 ln 

2B§m 2 (m - m s ) 4 M 2 2ffgm 2 (m - m a ) 4 M 2 - 2ffg (m - m a ) 4 (m + 8m s ) 2 M 2 
+ 729vr 2 M8 n 7^ + 2187tt 2 M8 n "T 2 " + 19683vr 2 M8 R 7^" 

2048.Bg (m - m s ) 4 (m + 8m s ) (2m + m s ) 
+ 6561M| ^ 

20485^ (m - m s ) 4 (m + 8m s ) (m + 2m s ) 
+ 19683M8 

10245^ (m - m s ) 4 (l6m 2 + 109mm s + 64m, 2 ) r 

6561M| Le(M) 

40965^ (m - m s ) 4 (4m 2 + 7mm s + 16m 2 ) 

6561M8 Ll 



40965 2 (m - m s ) 4 (2m 2 + 9mm s + 16m 2 ) 2B% (m - m s ) 4 (37m 2 + 16mm s + 64m 2 ) 
6561M8 Ls + 19683vr 2 M8 

Diagrams 6 

^m 2 ln^ + ^(m + m s ) 2 ln^ + -^(m + 2m s ) 2 ln^ (43 ) 

727T 2 fl 2 2167T 2 /i 2 19447T 2 M 

+ 1944^ (3W + 22mms + 13m «) 
Diagrams 7 

£ 3 m 2 (m - m a ) M 2 i? 3 m, (m 2 - m 2 ) M 2 

27tt 2 M 2 ii 2 81vr 2 M 2 /i 2 1 J 

(m — 4m s ) (m — m s ) (m + 2m s ) M 2 B 2 (m — m s ) 2 (17m s + 28m) 

729M 2 tt 2 n 7^ 729vr 2 M 2 

Diagrams 8 

2i? 4 m 2 (m — m s ) 2 M 2 2B 2 m 2 (m - m a ) 2 M 2 K 

81vr 2 M 4 fi 2 ^ 243vr 2 M 4 ^ 2 { ' 

2B 2 (m-4m s ) 2 (m-m s ) 2 M 2 2B 2 (m - m s ) 2 (28m 2 - 8mm s + 25m 2 ) 
+ 2187vr 2 M 4 2187tt 2 M 4 
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• Diagrams 9 



Bq-iti 2 (m — m s ) 2 M 2 Blm(m — m s ) (m + m s ) M\ 

81vr 2 M4 ~iJ 243vr 2 M4 ^ 46 ^ 

o 

Bq (m + 8m s ) (m — m s ) 2 (m + 2m s ) M 2 Bq (m — m s ) 2 (m + m s ) 2 Mj^ 
+ 2187tt 2 M4 n 7^" 162vr 2 M4 n 7^" 

i?0 ( m — w s ) 2 (l9m 2 + mm s + 16m 2 ) 



+ 



Diagrams 10 



2187vr 2 M4 



ABqIji 2 (m — m s ) 3 M 2 4Bgmm s (m — m s ) 3 

243tt 2 M6 T 2 " 729tt 2 M6 ( 7) 

4i?g (m — 4m s ) (m + 8m s ) (m — m s ) 3 

6561vr 2 M6 n 

O 

2B 5 (m + m s )m s (m-m s ) 3 M\ ±B\ (m - m s ) 3 (28m 2 + 13mm s - 32m 2 ) 
243vr 2 M6 n 6561vr 2 M6 

Let us denote C(m,m s ) the sum of all these diagrams. One could naively think that C is the 
NLO expression for C4 once the propagators are reexpressed in terms of the physical 77 mass, so that 
we just have to introduce a HO remainder dc° le of the form 

c pole = (J + c Pole d pole _ (4g) 

However, like in the case of the topological susceptibility \ [H]> one notices that C does not vanish in 
the limit m — > and/or m s — > 0, contrary to the expectation that c^ ole then, since the effective 
potential becomes independent of 6 in these limits. This definition of dc° le must thus exhibit a 
divergent contribution in these chiral limits, so that C4° le dc4 le can compensate the non-vanishing 
value of C to fulfill cf 16 -> 0. 

Such a behaviour of dc° le is not very satisfying, as we would like HO remainders to become smaller 
and smaller in the limits of vanishing quark masses. This can be resolved by defining the remainder 
as 

cf lc (m, m s ) = C(m, m s ) - C(m, 0) - C(0, m.) + c p 4 ole d^ le . (49) 

This corresponds formally to a redefinition of the HO remainder, but has the important consequence 
that dc° lc has no divergences and should remain small when approaching in the two massless limits. 
This remainder can thus be approximated by a polynomial, and in the physical setting where m s is 
much larger than m, we expect dc°^ e = 0(m 2 ). 

The difference between the two expressions Eqs. (j33|) and ()49[) corresponds to a reshuffling between 
the sum of LO and NLO contributions and HO remainders. Since we hope to achieve a formulation 
with small HO remainders, the numerical differences between these expressions may become impor- 
tant to extract information on the pattern of N = 3 chiral symmetry breaking from precise lattice 
values of C4 |j. For the moment, since no such precise determinations of C4 are available, we will not 
pursue the comparison further, and we will restrict our discussion to Eq. (|32p . bearing in mind that 
similar comments should certainly hold for the other formulations Eqs. (|33p and (|49p . 



6 Let us mention that in the simpler case of the topological susceptibility, the two formulations yield very comparable 
results in practice, as shown in Ref. [5]. 
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4 The cumulants of the topological charge and the pattern of chiral 
symmetry breaking 

4.1 Expected sensitivity to the three-flavour chiral condensate 

We are now in a position to discuss the sensitivity of the cumulants of the topological charge to the 
pattern of chiral symmetry breaking in the limit of three massless quarks. We see that the topological 
susceptibility Eq. ([27|) involves the three-flavour quark condensate 

E(3) = - lim (0|mi|0) = F$B , (50) 

m u ,m d ,m s ^0 

together with the Zweig-rule violating NLO LEC Lq in the combination [9j 

X o F$B + 32BiL r 6 s , (51) 

where s = m u + + m s . As discussed in Refs. |22j . for physical quark masses, this combination is 
actually very close to the quark condensate defined in the two-flavour chiral limit m u ,md 0, but 
m s kept at its physical value 



£2 ) = - lim OHO =^o + 32m s B^ M __^l g^°-^-^log— + 52 



where the ellipses denote HO contributions starting at O(m^). For simulations close to the physical 
case where m is much smaller than m s , the topological susceptibility essentially probes the two- 
flavour quark condensate (as the two quantities differ only by a contribution suppressed by a factor 
m/m s ), and consequently, x does not not probe the three-flavour quark condensate (defined in the 
corresponding chiral limit m u ,m c i,'m s — > 0), contrary to what is often asserted. Indeed, these two 
condensates can be significantly different in the case of large contributions from sea-quark ss pairs 
(encoded in the Zweig-rule violating LEC Lq), as first discussed in Ref. [22], and supported by 
several analyses (ttK scattering [23, 24J , scalar form factors [251428] . lattice data [9"ll29]). Therefore, a 
determination of "the" quark condensate from the topological susceptibility for lattice quark masses 
close to the physical ones can only be interpreted as a confirmation of our current knowledge on 
the pattern of N = 2 chiral symmetry breaking probed experimentally through irir (re)scattering 
and through lattice simulations of the spectrum of pseudoscalar mesons. But this does not shed 
much light on the pattern of N = 3 chiral symmetry breaking, and in particular on the size of the 
three-flavour condensate. 

This feature can be seen most easily once the topological susceptibility is expressed in the frame- 
work of Resummed Chiral Perturbation Theory (Re%PT) p [22l[M[27H29l [i3tii5] . Such a framework 
has been designed to cope with the fact that three flavor chiral series are not necessarily saturated 
by their leading-order term. Indeed if there are large contributions from sea-quark ss pairs and a 
suppression of the three-flavour condensate, a numerical competition between leading and next-to- 
leading orders will occur. In this case an arbitrary function of pseudoscalar observables will not 
necessarily have a well-convergent series and one should not truncate chiral series by neglecting 
higher-order contributions. 

In that framework the topological susceptibility for physical quark masses is given by 



no pole 

X ' 



FlMl r 



physical 2 2r + 1 



(7r + 2)e(r) 9 r SFjMj - \F l K M l K + F^M. 

2r(2r + l) + 2 (r - l) 2 (r + 2) F%M% 



+... (53) 



where e(r) is a simple function of r = m s /m defined in Eq. (|86|) . and the ellipsis denotes higher-order 
remainders starting at next-to-next-to-leading order (see Ref. [9] for a more extensive discussion of 
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the topological susceptibility) 0. It is clear from this expression that x is essentially sensitive to r, but 
exhibits no sensitivity to S(3). This is related to the fact, shown in Ref. [22], that the three- flavour 
chiral expansion of F 2 M 2 yields a strong correlation between the two-flavour condensate and r. 

If determining the three-flavour condensate from the topological susceptibility around the physical 
point seems hopeless, the situation becomes more promising when turning to the fourth cumulant 
C4, which involves Lq in two different combinations 

_ [31 

04^^50+3250^, -965^^-, (54) 

For simulations close to the physical quark mass hierarchy, the first term corresponds to the two- 
flavour condensate (like for the topological susceptibility), but the second term opens the possibility 
to disentangle Lg and S(3) by considering the dependence of C4 on the light-quark mass m. Therefore, 
contrary to the topological susceptibility x, the fourth cumulant C4 may provide a probe of the pattern 
of three-flavour chiral symmetry breaking even for simulations close to the physical case. 

Before turning to the numerical analysis of lattice settings close to the physical quark mass 
hierarchy, let us emphasise that simulations with degenerate quarks could be much more interesting 
to extract the three-flavour condensate from topological quantities. Indeed, in the limit where the 
iV quarks in Eqs (|27|) and (|31|) are degenerate with mass rh (leading to pseudoscalar mesons with a 
common physical mass M), the following combination 

v no pole 1 ^ 2 no pole, 3m 2 F 2 N B N N 2 - 1 3m 2 B% N 2 -l3rh*B% 2rhB N 3 

does not involve NLO LECs, and is particularly suited to determine the quark condensate S(iV) = 
F^Bj\f in the chiral limit of N massless flavours from simulations with N degenerate light quarks. 

4.2 Numerical analysis 

Several lattice computations of the (normalised) fourth cumulant have been performed in pure gauge 
theory for different number of colours [301432] . Even though their values for 3 colours are rather close 
to those that we will obtain in the following, they cannot be used to investigate the role of strange 
sea-quark pairs in chiral symmetry breaking (see Ref. [8] for more details and references on these 
simulations). There exist also computations with 2 and (2+1) dynamical fermions from the TWQCD 
collaboration (33HM], whose central values are in the same ball park as our results, but with so large 
uncertainties that any practical comparison with our figures is meaningless at the current stage. 

Our analysis will thus be limited to a proof of principle that the fourth cumulant C4 is indeed an 
interesting probe of the pattern of three-flavour chiral symmetry breaking. We consider Eq. (|33|) in 
the case of lattice simulations with (2+1) dynamical quarks of masses rh and rh s measured as 

m s rh 
P = — , Q = — ■ (56) 
m s m s 

We introduce the parameters involved in the LO chiral Lagrangian 

*(3) = m ZV-m. (57) 



F 2 M 2 ' W p2 



m 



One can also introduce the ratio of order parameters 



_ *(3) _ 2mB , , 

Y{3) ~Mf ■ (58) 



7 We take this opportunity to point out that Eq. (26) in Ref. [9] should be written with a factor 16 (rather than 4) 
in front of the term with L\ and Lg, so that the last two terms in this equation cancel and yield Eq. ()53|) , 
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Figure 4: C4 (in units of 10 4 GeV 4 ) as a function of X{3) and Y(3) for p = l,r = 1/q = 25 ; and aZZ 770 
remainders set to zero. 



X(3), Y(3) and Z(3) indicate the saturation of the three-flavour chiral expansions of F%M%, M 2 and 
F 2 by their leading-order contribution respectively. As indicated before, there are several indications 
that -X"(3) and Z{3) are not particularly close to 1, and that the three-flavour chiral expansions of 
pseudoscalar masses and decay constants experience a numerical competition between leading and 
next-to-leading orders. 

Denoting X an observable measured on the lattice, and X the same observable at the physical 
point, we obtain for the fourth cumulant on the lattice 



-no pole = F 2 M 2 pg (g 3 + 2) rX(3) M*p 2 q 2 «13q 2 + 22g + 37)r 2 Y(3) 2 

° 4 2(g + 2) 4 + 96vr 2 ( (? + 2) 4 l ' 

16Lg(/z)Mfc 2 g (g 4 + 2q 3 + 6q 2 + 8q + l) r 2 Y(3) 2 _ 96(3L 7 + L r 8 (fi))M*p 2 q 2 (q 3 + 2) r 2 F(3) 2 
(9 + 2) 4 (q + 2f 

O 

9M 4 p 2 g 2 r 2 y(3) 2 M 2 3M*p 2 q 2 (4g 3 - 3g + 2) r 2 Y(3) 2 M 2 
+ 32vr 2 (g + 2) 4 + 32vr 2 (g + 2)5 

O 

M 4 3pV(9 + 1) VF(3) 2 . M 2 - M%p 2 q(q + l) 2 (<? 3 - 3g 2 + 3g + 2) r 2 Y(3) 2 M 2 



327r 2 (g + 2) 4 /i 2 ^ 32vr 2 ((7 + 2) 5 /i 



2 



M 4 pV(2g + l) 2 r 2 y(3) 2 1 M 2 | M*p 2 q (2g 4 - 8g 3 + 13g + 2) r 2 y(3) 2 M 2 pol rfe 

111 2 + 9S8^2/„ j_ 9^4 111 „2 +C 4 a c 4 



96vr 2 (q + 2) 4 /i 2 288vr 2 (g + 2) 4 fi 

In Resummed xPT |24 [ l29 t l43]. Lg 78 (/x) are replaced by their expressions derived from the chiral 
expansion of observables expected to have good convergence properties [Fp and FpMp with P = 
n,K,rj). The corresponding formulae were given in refs [£29j and are recalled in App(Bj d Ci is the 
HO remainder which is set to zero in the present analysis (it is expected to be of order 0(m 2 ) and 
thus with the scaling d C4 ~ p 2 d Ci for simulations where q = m/m s is small). 
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0.05 0.10 0.15 0.20 0.25 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

q P 

Figure 5: C4 as a function of q — rh/rh s for p = 1 (left panel) and as a function of p = fh s /m s for q = 0.15 
(right panel). Light (green), dark (blue) and black curves correspond to X(3) = 0.2, 0.5, 0.8 respectively, whereas 
solid and dashed curves correspond to Y(3) = 0.9 and 0.2 respectively. We take r = 25 and all HO remainders 
are set to zero. 




Figure 6: x fl e ft panel) and 62 = Ci/(12x) (right panel) as a function of q = m/m s for p = l,r = 25. Light 
(green), dark (blue) and black curves correspond to X(3) = 0.2,0.5,0.8 respectively, whereas solid and dashed 
curves correspond to Y(3) = 0.9 and 0.2 respectively. All HO remainders are set to zero. 



The lattice counterpart of Eq. (|32H can be obtained by replacing all masses in the logarithms by 

o 

their LO values denoted as Mj>, yielding 

-no pole _ F 2 M 2 pq (g 3 + 2) rX{3) M*p 2 q 2 ((13q 2 + 22q + 37)r 2 F(3) 2 

° 4 2((? + 2) 4 + 96vr 2 (g + 2) 4 1 j 

mU G {ri)M^p 2 q (q A + 2q 3 + 6q 2 + 8q + l) r 2 y(3) 2 96(3L 7 + L r 8 (n))M*p 2 q 2 (q 3 + 2) r 2 Y{3) 2 



(q + 2f (q + 2f 

O O 

3M*pV (2 + q 2 ) r 2 Y{3) 2 M% M^p 2 q{q + l) 2 (q 3 +9q + 2) r 2 Y(3) 2 M 2 K 

+ 8vr 2 (g + 2)5 32vr 2 (g + 2)5 ~JT 



2 



Myg (g 4 + 2q 3 + 6q 2 + 8q + l) r 2 Y{3) 2 M*_ „ ao pole - pole 
1447r 2 (g + 2) 4 ^ 2 + ° 4 



We will use Eq. (|6Up to determine the potentiality of the fourth cumulant C4 to extract the 
three-flavour quark condensate from simulations performed with a mass hierarchy similar to the 
physical situation (for this first numerical investigation, we use = 139.6 MeV, Mr- = 495.7 MeV, 
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0.0 
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X(3) 



0.8 
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Figure 7: Order-by- order contributions to x (left panel) and C4 (right panel) as a function o/X(3) in the case 
p = l,q = 0.15, r = 25. Light (green), dark (blue) and black curves correspond to LO contributions, NLO ones 
and their sum respectively. Solid and dashed curves correspond to Y(3) = 0.9 and 0.2 respectively (they are 
identical for the LO contribution which depends on X(3) and quark masses only). All HO remainders are set 
to zero. 



M v = 
zero. 



547.8 MeV, F v = 92.2 MeV, F K /F n = 1.18). All higher-order (HO) remainders are set to 



We first consider the physical case where p = 1, q = 1/r, and in Fig. HJ we show the dependence of 
C4 on X(3) and Y(3) assuming r = 25. We see a noticeable sensitivity to X{3) and a weak sensitivity 
to Y(3), as expected from the discussion in Sec. [H Let us move from the physical case to consider 
simulations with a near physical strange quark mass and light quark masses larger than the physical 
one (i.e., p close to 1 and q > 1/r). The variation of C4 with respect to p and q, X(3) and Y(3) is 
shown in Fig. The slope of the curve is very dependent on X(3), and more weakly on Y(3), and 
the dependence of C4 on the light-quark mass rh provides an efficient probe of the three-flavour quark 
condensate X(3). 

As indicated in the qualitative discussion of the previous section, the topological susceptibility \ 
has a much weaker sensitivity to -X"(3), which arises only for values of q closer to 1, i.e., light quark 
masses away from the physical situation. This can be seen from the left-hand side of Fig. [U] (we 
plot here Eq. (f28|) . obtained without singling out the r/-pole and denoted x n ° pole i n Ref. [9])- We 
also display the normalised fourth cumulant 62 = C4/(12%) on the right hand-side of Fig. which 
exhibits the same appealing sensitivity to the three- flavour quark condensate as C4. We notice that 
its dependence on the simulated light quark mass is modified compared to the LO result [35J 



52;LO 



2 + q A 



12(2 + qf 



(61) 



both for the value at the origin and the slope when the quark condensate does not saturate the thee- 
flavour expansions of quark masses. Different patterns of three-flavour chiral symmetry breaking, 
and in particular, different values of the three-flavour quark condensate, can modify the dependence 
of this cumulant on the simulated light quark mass in a striking way. 

In Fig. we indicate the LO and NLO contributions to x an d C4 for a simulation with p = 1, q = 
0.15 and r = 1/q. A significant cancellation of the two contributions in the case of C4, emphasising the 
need to perform the analysis of this quantity in xPT up to one loop, whatever the pattern of chiral 
symmetry breaking. Moreover, this partial cancellation does not prevent the sum of LO and NLO 
contributions for C4 to exhibit a strong dependence on X(3), contrary to the topological susceptibility 
X which again shows a very limited sensitivity to X(3) (as expected since q = rh/fh s is small). 

Finally, we illustrate typical values for the fourth cumulant C4 and 62 for current lattice simula- 
tions. We consider the results of fit B5 from Ref. [3], corresponding to the analysis of the spectrum 
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r 


25.6 ± 1.0 


X(3) 


0.55 ±0.04 


Y(3) 


0.91 ±0.07 


Z(3) 


0.60 ± 0.04 




1.25 ±0.02 


m s (2 GeV)[MeV] 


81.2 ±6.0 


m(2 GeV)[MeV] 


3.2 ±0.3 


So /3 (2 MeV)[GeV] 


243 ± 12 


S (2 GeV) [GeV] 


2.81 ±0.38 


FolMeV] 


71.3 ±2.4 


L 6 (M p ) • 10 3 


0.04 ±0.07 


L 7 (M p ) • 10 3 


-0.25 ±0.14 


L 8 (M p ) • 10 3 


1.04 ±0.09 



Table 1: Results taken from Ref. [91 obtained through a fit to the RBC/UKQCD data on pseudo scalar masses, 
decay constants and topological susceptibility JJEjj within the framework of Resummed Chiral Perturbation 
Theory. 



Case 


P Q 


X • 10 4 [GeV 4 ] 


c 4 ■ 10 4 [GeV 4 ] 


b 2 


L s = 24 


1.71 


0.192 


5.66 ± 1.54 


-0.81 ±0.47 


-0.012 ±0.009 




1.71 


0.311 


8.98 ±2.61 


-0.87 ±0.75 


-0.008 ± 0.010 


L s = 32 


1.61 


0.154 


3.88 ±0.83 


-0.67 ±0.32 


-0.014 ± 0.007 




1.61 


0.220 


5.60 ± 1.26 


-0.84 ±0.47 


-0.013 ± 0.008 




1.61 


0.286 


7.23 ±1.71 


-0.89 ±0.61 


-0.010 ± 0.009 


Physical 


1 


l/r 


0.40 ± 0.02 


-0.10 ±0.01 


-0.021 ± 0.004 



Table 2: Predictions for the topological susceptibility and the fourth cumulant (without identifying the r\ pole) 
for different lattice data sets produced by the RBC and UKQCD collaborations \4@1 and in the physical case, 
following the analysis performed in Ref. leading to the results summarised in Table 1. 



of light pseudoscalar mesons and the topological susceptibility from the RBC and UKQCD collab- 
orations [3B] within the framework of Resummed xPT [2411291135] E|. This analysis led to the values 
recalled in Table [TJ where one notices the rather low values for the quark condensate and the pseu- 
doscalar decay constant in the three-flavour chiral limit: X(3) = 0.55 ± 0.04 and Z{3) = 0.60 ± 0.04. 
Using the outcome of this fit (including the correlations among the LO order parameters and the 
HO remainders), we obtain the predictions for Xi c 4 an d ^2 collected in Table [2] for the different sets 
considered in Ref. [9j as well as in the physical case. The negative value of 62 indicates that the 
distribution of the topological charge is flatter than a pure Gaussian. The values for p = rh s /m s 
are different from those quoted in Ref. [36], since we reassessed the determination of the lattice 
spacing and quark mass based on the mass of the O baryon in ref. [9]. Finite- volume effects are not 
included, whereas the effect of HO remainders is taken into account through the fit and an additional 
uncertainty coming from d Ci ~ 0.15 has been added to the results for £4 and 62- As in the case 
of the topological susceptibility illustrated in Ref. [S], the impact of d CA can be easily modeled and 
determined from a fit once a sufficiently large set of values at different simulated quark masses is 
available. Obviously, it would be very interesting to compare these predictions for C4 and 62 with 
data from lattice simulations. 

8 A recent update of the lattice results has been made by the RBC and UKQCD collaborations in ref. [37], confirming 
our analysis concerning the difficulties of a precise determination of lattice spacings and adding a new set of data, which 
is not included to produce the predictions shown here. 
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5 The two-flavour case 



As indicated after Eq. (|3ip . it is possible to compute the topological susceptibility or the fourth 
cumulant using the N = 2 xPT Lagragian described in Ref. [lj and built around the two-flavour 
chiral limit m u = = (but m s kept physical). The only degrees of freedom of the theory are 
soft pions and the resulting expansions are series in powers of m u and nid only. The LECs l\ and hi 
involved are by construction different from those in N = 3 xPT (denoted Li and Hi), as they are 
defined in two distinct chiral limits. For instance, at leading order, two-flavour xPT will involve the 
pseudoscalar decay constant F and the quark condensate —F 2 B defined in the limit m u = nrid = 0, 
but m s physical, and one has 

lim F 2 = F 2 , lim -F 2 B = -F$B , (62) 

showing that the two sets of quantities are in principle different: one expects that two-flavour chiral 
order parameters are larger than the three- flavour ones [22] , and our fits to lattice data indeed confirm 
this trend with a significant suppression of chiral order parameters when the strange quark mass is 
sent to zero pi29]. 

A further complication stems from the fact that the structure of NLO chiral Lagrangian used in 
Refs. [2] and pQ is different, so that the correspondence between the LECs in the two theories is not 
immediate beyond leading order. It is straightforward to repeat the same arguments as in Sec. 12.21 
with the generating functional of N = 2 xPT given in Ref. [1], in order to obtain the two-flavour 
expressions in the isospin limit 

X nop ° le = ^^ + 2 m 2 B 2 [i^) + h^)-£ 7 -h 3 ]- 3 -^lo g ^ + 0(m 3 ) (63) 

cT P ° le = -^^+ 9 -^-2rn 2 B 2 [il^) + h^)-i 7 -h 3 } (64) 



3m 2 B 2 

' 1287T 2 



, 2mB ol M 2 
log — 2~ + 31og^- 



+ 0(m 3 



where we have performed the same separation between unitarity and tadpole contributions for C4 
as in the iV = 3 case. These formulae should be used to extract the two-flavour chiral condensate 
from simulations performed at m s = m s and various values of the light quark mass rh. One notices 
in particular that the sum of Eqs. ([63]) - ([65]) . i.e. C4 + x, yields a result in agreement with Eq. ([55]) : 
the cancellation of all NLO LECs shows that this sum would be a particularly clean probe of the 
two-flavour condensate —BF 2 in the isospin limit. 

One can go one step further and actually match N = 2 and N = 3 chiral theories in order to 
exhibit the m s -dependence of the N = 2 LECs [1,48-50], yielding the NLO matching formulae for 
the quark condensates and pseudoscalar decay constants 

F 2 = F 2 + 16 ms B Ll(v)-^log^ (65) 

F 2 B = F 2 B + 32m s B 2 Ll^)- 1 ^log^-^-log^^, (66) 

as well as for the relevant NLO LECs 



2m 2 B 2 [e 3 ^) + h[(p)-£ 7 -h 3 ] = 24m 2 J B 2 [^ + 3^ 7 + ^]-^log^ (67) 



m 2 B 2 4m s B 5m 2 B 2 F 2 B m 2 



log 



96vr 2 b 3/x 2 144vr 2 4m s 
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where -^7 8 = ^6 7 8-N=3 denote the LECs from the N = 3 chiral Lagrangian. It is easy to check 
that Eqs. ()63[) - (|65() correspond indeed to the second-order expansion in m of the three-flavour results 
Eqs. ([28]) and (|33|) . providing a further cross-check of these equations. 

At first sight, one might be surprised that two high-energy counterterms h\ and I13 arise in these 
expansions. Indeed, these high-energy counterterms encode physics of higher energy scales that are 
not included dynamically in the theory and their value is dependent on its ultraviolet regularisation. 
This seems to contradict the fact that x an d C4 are topological quantities related to chiral symmetry 
breaking arising in QCD at low energies. This paradox can be solved thanks to the matching with 
three-flavour xPT. The two-flavour chiral expansions of x an d C4 involve only the difference of high 
energy counterterms h\ — hz, which can be matched with three- flavour xPT as in Eq. (|67p . The NLO 
expansion of h\ — /13 in powers of m s involves only N = 3 low-energy constants, but none of the high- 
energy counterterms Hi. It means that h\ — /13 is a combination of N = 2 high-energy counterterms 
which is characterised by the dynamics of K and rj mesons but is independent of the ultraviolet 
regularisation of the theory in relation with more massive (non-Goldstone) degrees of freedom. Indeed 
we have seen that our three-flavour results for x an d C4 rely crucially on the propagation of the 77- 
meson, which is not a dynamical degree of freedom of two-flavour xPT- Therefore, all the diagrams 
in N = 3 involving the propagation of K and 77 mesons (such as Fig. [1]) must be absorbed into 
high-energy counterterms once the computation is performed in N = 2 xPT. The presence of h\ — /13 
in Eq. (|65p is thus normal and can be easily explained by the peculiar role of the rj meson in the 
NLO expansions of x an d C4. 

Finally, let us stress the ambiguity of the alternative expression used in the literature [35 , 46J and 
derived from Eqs. (|27p and (|3ip by setting N = 2. The resulting expressions have the same structure 
as Eqs. ([63]) and ([65]) . but they involve LECs labeled Lq, Lj, L& that are very easy to misunderstand. 
Their names allude the three-flavour case, but they are actually defined in the two-flavour chiral 
limit. In our notation, we would denote them as £6,7,8;iV=2 which have definitions and values that 
differ from those in three-flavour %PT (which we would denote -^6,7,8^=3) • The connection between 
these LECs L6,7,8;iV=2 and the usual ti and hi can be obtained from refs. [UHUH2] 

£ r 3 = -8Ll. N=2 - 4L r 5 . N=2 + 16Lg ;iV=2 + 8L 8; JV=2 , (68) 

£ 7 = -16^=2-8^=2, (69) 

h\ = 8L1. N=2 + 4L r 5 . N=2 — 4L s -n=2 + 2H 2 -,n=2 , (70) 

hi = ALl. N=2 + 2H 2 . N=2 . (71) 

In order to avoid any confusion and make a direct link with phenomenological analyses performed in 
N = 2 xPT [1], one should always use the expressions Eqs. ()63p and (|65p to deal with the two-flavour 
chiral expansions of the topological susceptibility and the fourth cumulant. 

6 Conclusion 

Due to their connection with the ^Ta(I) anomaly, topological observables describing the distribution 
of the topological charge (or gluonic winding number) are able to probe the structure of low-energy 
QCD, and in particular the pattern of chiral symmetry breaking. Moreover, they can be accessed 
through lattice QCD simulations and are thus complementary to the analysis of the quark mass 
dependence of observables related to light pseudoscalar mesons. The two most prominent topological 
observables are the topological susceptibility x an d the fourth cumulant C4, corresponding to the 
first two non-trivial moments of the distribution of the topological charge. The present paper is a 
follow-up of the investigation started in Ref. [9j, in order to determine if these observables can be used 
to constrain the pattern of three-flavour chiral symmetry breaking (defined in the m u = rrid = rn s 
chiral limit). We have performed the computation of the two quantities up to next-to-leading order 
in three-flavour Chiral Perturbation Theory using two methods: the elegant approach of the effective 
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potential and the more pedestrian computation of Feynman diagrams. In the first case, we derived 
expressions valid for an arbitrary number of light flavours, whereas in the second case, we were able to 
pin down the role played by the propagation of flavour-singlet mesons. This allowed us to emphasise 
the importance of reexpressing the denominators arising at next-to-leading order expression and 
involving the leading-order contribution of the ?7-mass in terms of its physical value. 

We then investigated the sensitivity of the next-to-leading order expressions of the topological 
susceptibility and the fourth cumulant to the pattern of three-flavour chiral symmetry breaking in 
the context of lattice QCD simulations. We confirmed the result of Ref. [9] that for light-quark mass 
hierarchy close to the physical case, the topological susceptibility has only a sensitivity to the quark 
condensate in the two-flavour chiral limit (m u = rrid = but m s physical), but not to the three- 
flavour one. Since there are phenomenological indications that the two condensates are significantly 
different, one has to devise alternative ways of obtaining the three-flavour quark condensate from 
topological observables. We showed that the fourth cumulant has indeed a much better sensitivity to 
the three-flavour condensate than the topological susceptibility and was able to determine whether it 
is suppressed compared to the two-flavour one. As an illustration, we provided numerical predictions 
for the two topological observables in the case of the results of the RBC and UKQCD collaborations 
already analysed in Ref. [9]. 

In addition, we discussed the case of three degenerate quarks, and we showed that a combination 
of the topological susceptibility and the fourth cumulant is able to pin down the thee-flavour quark 
condensate in a particularly clean way, as all next-to-leading order contributions cancel in this com- 
bination. We also considered the formulae for the two topological observables in two-flavour Chiral 
Perturbation Theory, and discussed the convention and physical content of the low-energy constants 
arising in the resulting chiral expansions. 

From the viewpoint of effective field theory, our results show clearly that the fourth cumulant 
is particularly well suited to investigate the three-flavour quark condensate using lattice QCD sim- 
ulations, and that it should be preferred to the topological susceptibility which has only a limited 
sensitivity to this chiral order parameter for hierarchy of quark masses close to the physical case. 
Obviously, it remains to be seen if this quantity can be computed accurately using lattice techniques. 
But we hope that our results will prompt lattice collaborations to study the fourth cumulant of the 
distribution of the topological charge in more detail, providing a new handle on the chiral structure 
of QCD vacuum. 
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A Diagram contributions to C4 

As indicated in Sec. 13.11 there are several classes of diagrams contributing to the chiral expansion of 
C4 up to one loop. In addition to the five LO diagrams in Fig. Q] dressed by adding a counterterm 
or a loop at the level of the vertices or the propagators, there are five additional diagrams collected 
in Fig. [3l In the following, we will give the contribution of the various diagrams with the same 
distinction for the argument of logarithms coming from tadpole or from unitarity contributions as in 
Eq. PD- 
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Diagrams 1 
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Diagrams 2 
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o 
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o 

B 2 (m - m s ) (m 2 - llmm s - 8m 2 ) M 2 

1458vr 2 (2m s + m) 11 7?" 

1024Bg (m s + 2m) (m - m s ) 2 , 256B 2 (8m 2 + 29mm s + 8m 2 ) (m - m s ) 2 

-L 6 [(i) H — — L 7 



81(m + 2m s ) bUV 81 (m + 2m s ) 2 

256-B 2 (4m 2 + 15mm s + 8m 2 ) (m - m s ) 2 r 
+ 81 (m + 2m s ) 2 ^ 



Diagrams 3 
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o 
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216vr 2 (m + 2m s ) 2 n T* 2 " 

o 
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In 



972vr 2 (2m s + m) /x 2 
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2 
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Diagrams 4 
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m 
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Diagrams 5 
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• Diagrams 7 
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• Diagrams 8 

B 2 m 2 (m - m s ) 2 ^ M 2 £ 2 m 2 (m - m s ) 2 ^ £ 2 (m - 4m s ) 2 (m - m s ) 2 M 2 

187r 2 (m + 2m s ) 2 11 /i 2 + 547r 2 (m + 2m s ) 2 H /i 2 4867T 2 (m + 2m s ) 2 " /i 2 

Bq (in — m s ) 2 (28m 2 — 8mm s + 25m 2 ) 
+ 486vr 2 (m + 2m s ) 2 



Diagrams 9 
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• Diagrams 10 
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Let us notice that vertices from L\ with derivatives applied to the internal propagators generate 
contact terms which yield tadpole contributions even in the case of the scattering diagrams 9 and 
10 where such contributions would not be expected naively. One can easily obtain the equivalent 

o 2 

contributions with the prescription of Eq. (|32p by setting the LO masses Mp in the argument of all 
logarithms. 



B Low-energy constants in Resummed Chiral Perturbation Theory 

In Resummed xPT [9"ll2"4" l l2T? | l4"5] ■ the chiral expansions of Fp and FpMp are expected to have small 
higher-order remainders when expressed in terms of the couplings arising in the chiral Lagrangian. 
These exact mass and decay constant identities can be inverted to reexpress NLO LECs in terms of 
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LO parameters of the chiral Lagrangian, physical quantities, and HO remainders 



F(3)AL 4 
Y(3)AL 5 
Y 2 (3)AL 6 
Y 2 (3)AL 8 
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so that e(r2) = 0, and e(ri) = 1 with r\ = 2{Fk Mk) / (F w M w ) — 1 ~ 8. d, d', e and e' are combina- 
tions of HO remainders associated with the chiral expansions of it, K masses and decay constants 
respectively. ALj = Ll(fi) — Li(fi) are independent of the renormalisation scale fj, and combine the 
(renormalized and quark-mass independent) constants £4,5,6,8 together with chiral logarithms: 
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Similarly, one can invert the relation for the r\ mass, yielding 



[y(3)] 2 L 7 



F 2 



32(r - l) 2 Ml 
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(91) 



where F^ is not measured, but can be computed using 



F 2 = F 2 Z(3)+8(r + 2)Y(3)M 2 AL 4 + -(2r + l)Y(3)M 2 AL 5 

M x M u 



(92) 



+ 



and we have introduced the difference of HO remainders 
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